Critical magnetic fluctuations induced superconductivity and residual density of 

states in CeRhln^ superconductor 
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We propose the multiband extension of the spin-fermion model to address the superconducting d- 
wave pairing due to magnetic interaction near critical point. We solve the unrestricted gap equation 
with a general d-wave symmetry gap and find that divergent magnetic correlation length £ leads 
to the very unharmonic shape of the gap function with shallow gap regions near nodes. These 
regions are extremely sensitive to disorder. Small impurity concentration induces substantial residual 
density of states. We argue that we can understand the large N res (0) = lhriT^o C V {T)/T value and 
its pressure dependence of the recently discovered CeRhIn$ superconductor under pressure within 
this approach. 

PACS numbers: 74.20, 74.20-z,74.50 



Recent discovery of,.-, superconductivity in 
CeMIn 5 (M = Co, Rh, Ir)[M has spurred renewed in- 
terest of heavy fermion systems and about the nature of 
its superconductivity. The rich phase diagram of these 
compounds and the tunability by pressure ancLchemical 
substitution of the transition metal elements] provide 
valuable information about the compctction/interplay 
between the magnetism and superconductivity and per- 
haps a possible quantum criticality as a unifying origin 
of the phase diagram of these compoundsQ While the 
complete understanding of the phase diagram and the 
underlying mechanism for the magnetism and the super- 
conductivity is still lacking, there are many details of 
thermodynamic and transport properties of each phase 
in these materials which need to be understood: large 
remnant iV res (0) — lmvr_o C P {T)/T values both in mag- 
netic and in superconducting phasesp, the weak [first 
order transition to superconductivity in CeRhIn$l, a 
unidentified new phase inside the mixed stateO, a strong 
deviation of ACiXf)/C(T c ) from BCS value and its pres- 
sure dependenceO'Q etc. — 

Motivated by CeRhIn$ experiment^, in this paper we 
examined possible conditions for a d-wave superconduc- 
tors (1) to create a substantial density of states (DOS) in 
superconducting phase with a small impurity concentra- 
tion, and (2) to have a large variation of it as a function 
of pressure while keeping the constant T C U. 

The N res (0) value at lQ.bkbar seeni-experimentally is 
almost half oi the normal state AT(0)B. From the pre- 
vious studiesu, it is required to have a large amount 
of impurities in order to create large iV res (0) by mag- 
netic/nonmagnetic isotropic impurity in the supercon- 
ducting state with lines of nodosa (the normal state scat- 
tering rate V should be about T/Ao > 0.5 with Born scat- 
terer and T/Ao > 0.1 with a unitary scatterenl3). Even 
more intriguing is the pressure dependence of N res (0); it 
varies from almost the half of the normal state N(0) at 
16.5 kbar to almost zero value at 21 kbar in one sample. 



Assuming a simple form of d-wave paring Ao(cos(fc x ) — 
cos(fcy)), such a large variation of N res (0) in the same 
sample requires extremely sharp increase of Ao with in- 
creasing pressure. This would be difficult to reconcile 
with almost constant T c values between 16 ~ 21 kbaru. 

To address these questions we need a microscopic 
model for the superconducting state in CeMIn^ materi- 
als. From the phase diagram of these materials 1-5 and 
from the thermal conductivity measurement indicating 
the unconventional pairing with lines of nodesO, we ar- 
gue that it is highly plausible that the superconductivity 
is mediated by the magnetic fluctuationaMLH. We propose 
a multiband generalization of the spin-fermion modeO 
where localized Ce spins S are interacting with the con- 
duction electrons (predominantly d band of In) via the 
Kondo exchange coupling J. In mixed momentum and 
real space representation the Hamiltonian is written as 

H = 4 (k)e(k)c Q (k) + jg (r) • 4(r)? a/ , C/ ,(r) + H s 

(i) 

where the first term is a kinetic energy and the sec- 
ond describes the Kondo exchange between Ce spins and 
conduction electron spin density. The last term repre- 
sents an effective low energy Hamiltonian for the local- 
ized spins. The dynamics of the localized spins without 
long range AFM order coupled with conduction elcGtrons 
is well captured by the spin correlation functioriljii 5 ], 



X(q^) = MS»(q, W )S'(q,-u,)) = i»,„ 0+ A-Qr+i > 
where loq is a spin relaxation energy scale, Q is the 2- 
dimensional antiferromagnetic vector, and £ is the mag- 
netic correlation length. The physics of this modcliop-a 
one band case has been investigated for a long timdlirtL 5 ! 
It is important to mention that the spin-fermion model of 
our multiband case is different in the following aspects: 
(1) the effective spin-fermion coupling should be much 
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weakerlll than that of the one band model. One impor 
tant consequence of it is that the relaxational energy, 
wo is much larger than that of the one band case w s /l 
(2) more importantly, while the £ 3> a limit can not be 
reached in the self-consistent one band model without 
deforming the FS topology of the conduction bandES, in 
the two band model in conjunction with the weak cou- 
pling above mentioned, the £ 3> a limit doesn't -neces- 
sarily modify the conduction band Fermi surfaced; (3) 
below T c , however, the spin fluctuation damping luq is 
suppressed for to < 2A SC , thus modifying the spin fluc- 
tuations spectrum. This effect should be dealt with in 
the fully self-consistent Eliashberg equation. In our weak 
coupling BCS approximation, this effect is ignored and 
justified by the fact A sc < ujq a posteriori. 

The phase diagram for this material! suggests that £ 
becomes gradually shorter away from the phase bound- 
ary (P c ~ 15kbar) toward the higher pressure. Then 
the basic feature of the magnetic fluctuations mediated 
potential is that near critical point with divergent £ 
V(q) = xil^ = 0) is sharply peaked around the an- 
tifcrromagnctic wave vector q ~ Q without much tail 
at other q values. With this potential the d-wave gap 
becomes mostly confined near the antinodal paints and 
the gap in the nodal regions becomes shallowEO. Upon 
decreasing the correlation length £, the potential has a 
longer tail extending outside the Q vector. As a result 
the gap function approaches a simple harmonic function 
~ (cos(k x )—cos(k y )) with a constant slope of the gap over 
most of the Fermi surface (FS). In this paper we show 
that the OP slope of the nodal region is very sensitive to 
the shape of the pairing potential (see Fig.l(a)-(b)) and 
the shallow gap near nodal points leads to the enhanced 
sensitivity to disorder so that only small amount of im- 
purities is required to produce substantial residual DOS 
N r ea (0) ■ To capture this detailed shape of the d-wave OP, 
we solved the unrestricted gap equation with only keep- 
ing the general d-wave symmetry. Then we study the im- 
purity effects using T-matrix approximation to calculate 
the impurity-induced DOS N res (uj) in the superconduct- 
ing state. 

Formalism. For CeRhln^ we assume Ao < uio- This 
allows us to use a weak coupling BCS gap equation of the 
spin-fermion model for the d-wave pairingtfl. For sim- 
plicity, we assume a circular FS in two dimensions and 
integrate out the perpendicular component of momen- 
tum up to the BCS cut-off energy ujd, which is naturally 
provided by o>o in our model. The effect of the impurity 
scattering is included within T-matrix approxiamtionS. 
For particle-hole symmetric case T3 = 0, and for d-wave 
OP with isotropic scattering Ti = (also without loss 
of generality we can choose T2 = by U(l) symmetry). 
Then we need to calculate only Tq(uj). Self energy is given 
E = rr , where T = uJttNq, N the normal DOS at 
the Fermi energy, n, the impurity concentration. Scatter- 
ing strength parameter c is related with the s-wave phase 



50 = ^0 Efc D2 +£ Y+A(fc) ' = W « + With this 

To the following gap equation is solved self-consistcntly. 
A(<f>) = -N J ^LVW -</>'). Ftf) 
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A(k) 



(2) 



Unlike the previous calculations of d-wave pairing, 
we assume no particular functional form for A(<^) ex- 
cept imposing D2 symmetry; namely A(wr/4) = 0(n = 
1, 3, 5, 7), A(4>) = A(<t> ± tt), and A(<f>) = -A(<f> ± tt/2). 
Therefore the gap equation can produce the most general 
d-wave symmetry gap solution for a given pairing poten- 
tial. The pairing potential is proportional to the static 
limit of x(q, w = 0) ~ ( g _Q^ + g-^ ; which is parameter- 
ized as follows. 



V{S<j>) = V d (b) 



((5^±7r/2) 2 + 6 2 ' 



(3) 



shift S as c = tan 1 (5). Now To(oj n ) 



30 On) 



where 



where the Lorenzian part is normalized and Vd(b) de- 
termines the total coupling strength. Apparently the 
parameter b is proportional to normalized in the 

circular Fermi surface (£ ~ an/b; a is the lattice param- 
eter). Both £ and Vd(b) are functions of pressure. We 
numerically determine Vd(b) tojpake T c constant in ac- 
cord with CeRhlns experiment^]. Finally, we introduce 
the FS weighting function F(<p) — cos^(2c/>) to correct 
the artifact of the circular FS aad to mimic the impor- 
tant aspect of real FS topologyL3. 

Results. In all calculations we use u>d = 1, which 
also serves the unit energy. In Fig. 1(a) we show the nor- 
malized pairing potentials V((f>)/Vd(b) as a function of b 
for illustration. In Fig. 1(b) the solutions of A(</>) for the 
Born limit scatterer (c = IT = 0.05) are shown for the 
potentials shown in Fig. I (a). 

The self-consistently determined Tb(cJ n ) is analytically, 
continued to real frequency using Pade approximant£3 
to calculate the self-energy (£o(w + i-v))- Then N(0) = 
-^2 k ImGo(cu,k) is calculated. In Fig.2.a-b we plot 
/mSo(w) for both Born and unitary limits for poten- 
tials with different b ~ £ -1 shown in Fig.l. There is 
no resonance at Fermi level for Born limit scatterer as 
known, but the scattering rate 7 = /mEo(w = 0) in- 
creases as the gap of nodal region becomes natter, and 
for all cases for \uj\ > Ao it approaches the normal scat- 
tering rate T N = T/(c 2 + 1) = 0.025. As a result the 
residual DOS N res (uj — 0) (see the Fig. 3. a) sharply in- 
creases with flatter gap region. On the other hand, for 
the case of unitary scatterer (Fig.2.b) there is a resonance 
at the Fermi level for all cases, but the strength of this 
scattering rate 7 has opposite trend in contrast to the 
Born limit; the flatter the gap is, the smaller 7 is. This 
is because the self-consistent equations are different for 
each case: 7 = Tuq/c 2 for Born limit and 7 = T/hq for 
unitary limit (no =< , 7 >). Also for \ui\ > A it 

^/7 2 +A(0) 
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approaches the normal scattering rate = T = 0.005. 
Nevertheless for both cases even small value of 7 is suf- 
ficient to create the substantial N res (0) at Fermi level 
when the gap is shallow near nodes. This is seen in 
Fig.3.a-b which plots the normalized DOS, N(u>)/Nq, for 
both Born and the unitary limits. The clearly seen trend 
is that the longer magnetic correlation £ produces a more 
residual DOS, N(lu)/No for a same amount of impurities. 

Fig. 4(a) summarizes this trend, shows plots of N(u> = 
0)/No as a function of b (~ P) for both impurity cases. 
The result shows that the Born limit scatterer has a 
stronger dependence on pressure compared to the uni- 
tary scatterer. This is because of the opposite trend of 
7 in unitary scatterer due to the resonant pole. In com- 
parison to the experimental data of N res (O)B(shown in 
Fig.4(b)), the Born scatterer fits the data better. 

In Fig. 5 we show T c suppression as a function of im- 
purity scattering parameter T for both Born (c = 1) and 
unitary (c = 0) limits for the representative potentials 
(b = 0.1). The impurity scattering parameters used in 
our calculations gives (T C (T) — T c o)/T c q suppression of 
at most a few % from Fig. 5. (unitary limit 0.8%, Born 
limit 4%). In passing we note r crit /A is about 50% 
larger than the simple cos(2<^>) type d-wave resultta both 
for Born and the unitary cases. This means that the un- 
restricted gap equation can find a more optimized gap 
solution up to higher temperature compared to the fixed 
form of A cos(2c6) solution. 

In summary, in this paper, we propose a multiband 
spin-fermion model description of the pairing in 

CeMIn^ materials. Motivated by the phase diagram of 
these materials, we assume that the magnetic correlation 
length £ decreases with pressure and the functional shape 
of the magnetic fluctuations mediated pairing potential is 
changing. Using this potential we show that the slope of 
the gap near nodes can be sharply changed. This strong 
change of the slope can explain the pressure dependence 
of N res (0)/No as well as the large value of it in CeRhln^ 
superconductorO close to the quantum critical limit of 
£ — > 00 with a small amount of impurity. 



FIG. 1. (a) The normalized pairing potential V((f>) /Vd(b) as 
a function of the exchange momentum cj>, for different values 
of 6(~ £ -1 )- In increasing order of potential width, 6=0.02, 
0.04, .... 0.2. (f> - it/2 is the AFM peak momentum Q. Inset 
is the Vd(b) which is numerically determined to make T c con- 
stant. The trend is that the potential height is decreasing as 
£ becomes shorter (or the pressure increases) and the width is 
increasing. (b)The OP solutions A((J>) for pairing potentials 
shown in (a) with impurities (Born limit c = 1 and T = 0.05) 
for all cases. In decreasing order of potential width, A(</>) 
becomes flatter near node. This trend of the flatter gap near 
node approaching the magnetic phase (the longer £) is also 
observed in the high T c superconductorscj. 
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FIG. 2. (a) The imaginary part of self-energy ImEo(w) for 
different pairing potentials shown in Fig.l. with Born limit 
scatterer (c = 1 and F = 0.05). With increasing potential 
width, 7 = MEo(w = 0) decreases.; (b) the same as (a) with 
the unitary scatterer (c = and T = 0.005). With increasing 
potential width, 7 increases. 
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FIG. 3. (a) The normalized DOS N{u)/N for different 
pairing potentials as shown in Fig. I. with Born limit scat- 
terer (c = 1 and F = 0.05). With increased potential width 
N(u))/N decreases. Inset: with /3 = 4 and T - 0.06. (b) the 
same as (a) with the unitary scatterer (c = and F = 0.005). 
With increased potential width N(lj)/No decreases. Inset: 
with P = 4 and F = 0.006. 
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FIG. 4. (a) The normalized DOS N{u = 0)/N as a func- 
tion of b (~ ~ P — P c ) for Born limit scatterer (c = 1, 
r = 0.05, solid square) and for the unitary scatterer (c = 0, 
T = 0.005, open square). Inset is with (3 = 4 (c = 1, T = 0.06, 
solid square; c = 0, V = 0.006, open square); (b) Experimen- 
tal data from Ref[5]. 
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FIG. 5. T c /T c0 as a function of r/A both for the Born 
(c = 1) and unitary limit (c = 0) with a typical paring poten- 
tial (b = 0.1). 

We are grateful to R. Movshovich, J. Thompson, J. 
Sarrao, Ar. Abanov for useful discussions. This work 
was supported by US DoE. Y.B. was partially sup- 
ported by the Korean Science and Engineering Founda- 
tion (KOSEF) through the Center for Strongly Corre- 



1 H. Hegger et al., Phys. Rev. Lett. 84, 4986 (2000). 

2 C. Petrovic et al., Europhys. Lett. 53, 354 (2001); C. Petro- 
vic et al., J. Phys. C ondens. Matter 13, L337 (2001). 

3 P.G. Pagliuso et al., |cond-mat/0107266 , 

4 N.D. Mathur et al., Nature 394 39 (1998); S.S. Saxena et 
al., ibid 406 587 (2 000). 



T7 



Rev. 
Phys. 



B 48, 4219 
Rev. B 50. 



5 R. A. Fisher et al., |sond-mat/010922l[ 

6 T. P. Murphy et al., |cond-mat/010417c 

7 E. Lengyel et al., unpublished. 

8 P. J. Hirschfeld and N. Goldenfeld, Phys 
(1993); R.Fehrenbacher and M. Norman 
3495 (1994). 

9 R. Movshovich et al., Phys. Rev. Lett. 86, 5152 (2001). 

If the data in Ref[5] were for a particular sample which hap- 
pens to contain a large amount of impurities, any cleaner 
sample should easily have twice as large T c . On the other 
hand the sample-to-sample variation of T c of CeRhlns is 
almost zero, so it means that the T c = 2.3K sample is a 
really intrinsic one. m 

1 The specific heat datatl suggests that the transition from 
the AFM to SC is weakly first order, we think that this 
weak first order transition doesn't change the nature of the 
AFM fluctuations drastically in the SC phase. 

2 Clear experimental evidence for the necessity of the two 
band model is that the resistivity of CeRhln^ continues to 
be metallic anrkeven decreased below Tjv in the magnetic 
phase (p < p c )o. It implies that the long range magnetic 
order of Ce spins doesn't induce the concomitant SDW gap 
in the conduction band. 

3 A.J. Millis, H. Monien, and D. Pines, Phys. Rev. B 42, 167 
(1990). 

4 P. Monthoux, A. V. Balatsky, and D. Pines, Phys. Rev. 
Lett. 67, 3448 (1991); Phys. Rev. B 46, 14803 (1992); K. 
Miyake, S. Schmitt-Rink, and C. M. Varma, Phys. Rev. B 
34, 6554 (1986). 

5 A. V. Chubukov, and D. Morr, Phys. Rep., 228, 355 (1997) 
and references therein. 

6 Jorg Schmal ian, D. Pines, and B. Stojkovic, cond 



mat/9804129 



Ar. Abanov, A. V. Chubukov, and A. M. Finkel'stein, Eu- 
rophys. Lett. 54, 488 (2001). 

Because of the circular FS, the pairing potential V(<f> — 4> ) 
itself is invariant with the simultaneous shift of 4> and cf> , 
and cannot distinguish the hot spots (parts of FS connected 
by the antiferromagnetic wave vector ±Q) and cold spots 
on the FS. Therefore it is essential to introduce the FS 
weight function to avoid the artifact of the circular FS. 
In this paper we use F(<p) = cos' 3 (20) with f3 = 8, but 
the results with /3 = 4 are not much different although it 
apparently produces less flat OP near nodes. 
H.J Vidberg and JW. Serene, J. of Low Temp. Phys. 29, 
179 (1977). 



4 



.-X.Shen (private communication) 



